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Abstract 


ShipMo3D  is  DRDC  Atlantic’s  object-oriented  library  for  modelling  ship  motions  in 
waves.  Previous  ShipMo3D  development  considered  ships  travelling  with  nominally 
steady  speed  and  heading.  This  report  describes  the  extension  of  the  ShipMo3D 
library  to  model  motions  of  freely  maneuvering  ships.  New  ShipMo3D  force  com¬ 
ponents  arise  from  hull  maneuvering,  resistance,  propulsion,  and  rudder-propeller 
interaction.  Comparisons  of  turning  circle  predictions  with  full-scale  trials  data  for 
the  tanker  Esso  Osaka  give  encouraging  results.  Comparisons  of  predictions  with 
motions  of  a  steered  warship  model  in  waves  give  very  good  results.  Excellent  agree¬ 
ment  between  predictions  for  a  freely  maneuvering  ship  and  for  a  ship  with  nominally 
steady  speed  and  heading  indicates  that  the  extension  of  ShipMo3D  to  freely  maneu¬ 
vering  ships  has  been  correctly  implemented.  It  is  recommended  that  future  work 
further  investigate  prediction  of  hull  maneuvering  forces. 


Resume 


ShipMo3D  est  la  bibliotheque  objet  utilisee  par  RDDC  Atlantique  pour  modeliser  les 
mouvements  des  navires  dans  la  houle.  La  version  precedente  de  ShipMo3D  se  limitait 
aux  navires  se  deplagant  avec  une  vitesse  et  un  cap  theoriquement  constants.  Dans 
ce  rapport,  nous  decrivons  les  developpements  de  la  bibliotheque  ShipMo3D  visant  la 
modelisation  des  mouvements  des  navires  manceuvrant  sans  contraintes.  Les  nouvelles 
composantes  de  ShipMo3D  simulent  la  force  exercee  sur  la  carene,  la  resistance,  la 
propulsion  et  les  interactions  gouvernail-helice.  Les  comparaisons  des  predictions  et 
des  donnees  tirees  d’essais  a  echellc  reellc  de  cercles  de  giration  par  lc  petrolier  Esso 
Osaka  sont  encourageantes,  tout  comme  les  comparaisons  avec  les  mouvements  d’une 
maquette  clirigee  de  batiment  de  guerre.  L’excellcnt  accord  entre  les  predictions  d’un 
navire  manceuvrant  sans  contrainte  et  un  navire  navigant  avec  une  vitesse  et  un 
cap  theoriquement  constants  indique  que  les  ameliorations  apportees  a  ShipMo3D 
visant  les  navires  manceuvrant  sans  contraintes  ont  ete  correctement  realisees.  Nous 
recommandons  que  les  travaux  a  venir  soient  consacres  a  une  etude  plus  poussee  de 
la  prediction  des  forces  exercees  sur  la  carene  lors  des  manoeuvres. 
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Executive  summary 

Simulation  of  Hydrodynamic  Forces  and  Motions  for  a 
Freely  Maneuvering  Ship  in  a  Seaway 

Kevin  McTaggart;  DRDC  Atlantic  TM  2005-071 ;  Defence  R&D  Canada  -  Atlantic; 
December  2005. 

Introduction:  In  support  of  DND  work  on  modelling  and  simulation,  DRDC  At¬ 
lantic  has  developed  the  ShipMo3D  library  for  simulation  of  ship  motions  in  waves. 
Previous  ShipMo3D  development  considered  ships  travelling  with  nominally  steady 
speed  and  heading  in  waves.  This  report  describes  the  extension  of  ShipMo3D  to 
prediction  of  motions  for  freely  maneuvering  ships  in  calm  water  and  in  waves. 

Principal  Results:  Forces  relevant  to  ship  maneuvering  have  been  introduced  to 
ShipMo3D.  Comparisons  of  turning  circles  with  full-scale  trials  data  for  the  tanker 
Esso  Osaka  give  encouraging  results.  Comparisons  of  predictions  with  motions  of  a 
steered  warship  model  in  waves  give  very  good  results.  Excellent  agreement  between 
predictions  for  a  freely  maneuvering  ship  and  for  a  ship  with  nominally  steady  speed 
and  heading  indicates  that  the  extension  of  ShipMo3D  to  freely  maneuvering  ships 
has  been  correctly  implemented. 

Significance  of  Results:  DND  now  has  an  object-oriented  library  for  modelling  of 
motions  of  freely  maneuvering  ships  in  calm  water  and  in  waves.  This  library  can  be 
used  to  develop  applications  in  support  of  acquisition,  operations,  and  training.  For 
simulation  of  motions  of  existing  ships,  agreement  between  predicted  and  observed 
motions  can  be  optimized  by  tuning  of  hull  maneuvering  force  coefficients. 

Future  Plans:  ShipMo3D  applications  and  corresponding  documentation  will  be 
developed  to  facilitate  simulation  of  maneuvering  ships  in  waves  using  input  ASCII 
files.  Documentation  will  also  be  developed  for  using  ShipMo3D  objects  within  simu¬ 
lations.  Future  work  will  also  likely  examine  improved  prediction  of  hull  maneuvering 
coefficients. 
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Simulation  of  Hydrodynamic  Forces  and  Motions  for  a 
Freely  Maneuvering  Ship  in  a  Seaway 

Kevin  McTaggart;  DRDC  Atlantic  TM  2005-071 ;  R  &  D  pour  la  defense  Canada  - 
Atlantique;  decembre  2005. 

Introduction:  En  appui  aux  travaux  de  modelisation  et  de  simulation  du  MDN, 
RDDC  Atlantique  a  cree  la  bibliotheque  objet  ShipMo3D  consacree  a  la  simulation 
des  mouvements  des  navires  dans  la  houle.  La  version  precedente  de  ShipMo3D  se 
limitait  aux  navires  se  deplagant  avec  line  vitesse  et  un  cap  theoriquement  constants 
par  rapport  a  la  houle.  Dans  ce  rapport,  nous  decrivons  les  developpements  de  la 
bibliotheque  ShipMo3D  visant  la  prediction  de  mouvements  des  navires  manoeuvrant 
sans  contraintes,  sur  une  mer  calme  ou  houleuse. 

Resultats  principaux:  Les  nouvclles  composantes  de  ShipMo3D  traitent  les  forces 
pertinentes  a  la  manoeuvre  des  navires.  Les  comparaisons  des  predictions  et  des 
donnees  tirees  d’essais  a  echellc  reellc  de  cercles  de  giration  par  le  petrolier  Esso 
Osaka  sont  encourageantes,  tout  comme  les  comparaisons  avec  les  mouvements  d’une 
maquette  dirigee  de  batiment  de  guerre.  L’excellcnt  accord  entre  les  predictions  d’un 
navire  manoeuvrant  sans  contrainte  et  un  navire  navigant  avec  une  vitesse  et  un  cap 
theoriquement  constants  indique  que  les  ameliorations  apportees  a  ShipMo3D  visant 
les  navires  manoeuvrant  sans  contraintes  ont  ete  correctement  realisees. 

Importance  des  resultats:  Le  MDN  dispose  maintenant  d’une  bibliotheque  objet 
pour  la  modelisation  des  navires  manoeuvrant  sans  contraintes  sur  une  mer  calme  ou 
houleuse.  On  pourra  utiliser  cette  bibliotheque  pour  elaborer  des  logiciels  en  appui  a 
L acquisition,  aux  operations  et  a  la  formation.  Lors  de  la  simulation  des  mouvements 
de  navires  reels,  on  optimisera  l’accord  entre  les  mouvements  predits  et  observes,  en 
ajustant  les  coefficients  des  forces  exercees  sur  la  carene  lors  des  manoeuvres. 

Travaux  ulterieurs  prevus:  Les  applications  de  ShipMo3D  et  la  documentation 
correspondante  seront  elaborees  afin  de  faciliter  la  simulation  de  navires  manoeuvrant 
dans  la  houle,  en  utilisant  des  hchiers  d’entree  en  ASCII.  On  redigera  egalement  la 
documentation  relative  a  l’utilisation  d’objets  de  la  bibliotheque  ShipMo3D  lors  des 
simulations. 
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1  Introduction 


Naval  modelling  and  simulation  applications  frequently  require  predictions  of  ship 
motions  in  calm  water  and  in  waves.  DRDC  Atlantic’s  ShipMo3D  is  an  object  oriented 
library  being  developed  to  meet  requirements  of  a  variety  of  simulation  applications. 
Previous  ShipMo3D  developments  [1,  2,  3,  4]  consider  the  forces  and  resulting  motions 
in  waves  for  a  ship  travelling  with  nominally  steady  speed  and  heading. 

The  present  work  extends  the  ShipMo3D  library  to  simulate  the  motions  of  a  freely 
maneuvering  ship  in  waves.  Force  components  arising  from  hull  lift,  hull  resistance, 
hull  cross-flow  drag,  rudder  lift  and  drag,  and  propeller  thrust  are  important  for  a 
freely  maneuvering  ship  and  are  incorporated  into  ShipMo3D.  Unlike  many  ship 
maneuvering  programs,  the  ShipMo3D  library  includes  retardation  forces  essential 
for  accurate  modelling  of  ship  motions  in  waves. 

The  next  section  gives  a  brief  literature  review  of  ship  maneuvering  prediction.  Sec¬ 
tion  3  gives  the  axis  systems  and  equations  of  motions  for  a  freely  maneuvering  ship, 
followed  by  discussion  of  hull  maneuvering  forces  in  Section  4.  Section  5  describes 
hull  cross-flow  drag  forces  and  their  relationship  to  the  hull  maneuvering  forces  in 
the  previous  section.  Hull  resistance  and  propulsion  are  described  in  Sections  6  and  7 
respectively.  Section  8  discusses  rudder  deflection  forces,  and  is  followed  by  treatment 
of  rudder-propeller  interaction  forces  in  Section  9.  Section  10  describes  modelling  of 
rudder  and  propeller  control  systems.  Section  11  gives  validation  of  turning  circle 
predictions  with  full-scale  trials  data  for  the  tanker  Esso  Osaka.  Modelling  of  maneu¬ 
vering  using  unclassified  data  for  the  Canadian  Navy’s  HALIFAX  class  is  described 
in  Section  12.  Section  13  gives  validation  for  predicted  motions  in  waves  using  ex¬ 
perimental  data  for  a  steered  warship  model.  Recommendations  for  future  work  are 
given  in  Section  14,  followed  by  final  conclusions  in  Section  15. 

2  Literature  Review  of  Ship  Maneuvering 
Force  Prediction 


Ship  maneuvering,  also  spelt  “manoeuvring”  and  “maneuvring” ,  has  been  the  subject 
of  extensive  literature.  For  application  within  ShipMo3D,  maneuvering  simulation 
approaches  of  greatest  interest  are  those  which  predict  ship  maneuvering  properties 
based  on  hydrodynamic  forces  acting  on  the  hull  and  appendages  (including  pro¬ 
peller),  and  aerodynamic  forces  acting  on  the  superstructure.  The  present  report  only 
considers  hydrodynamic  forces.  Incorporation  of  aerodynamic  forces  into  ShipMo3D 
will  likely  be  the  subject  of  future  work. 

Several  references  provide  good  overviews  of  ship  maneuvering  forces.  Bertram  [5] 
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presents  one  of  the  most  recent  comprehensive  discussions  of  maneuvering  forces.  His 
book  also  covers  other  relevant  areas  of  ship  hydrodynamics,  including  resistance, 
propulsion,  and  seakeeping.  The  International  Towing  Tank  Conference  Specialist 
Committee  on  Esso  Osaka  [6]  have  also  produced  a  good  recent  reference  on  ship 
maneuvering  forces,  with  emphasis  on  maneuvering  predictions  for  the  Esso  Osaka,  a 
tanker  that  has  been  the  subject  of  extensive  full-scale  and  model  test  studies.  Barr 

[7]  compares  several  different  programs  for  simulating  ship  maneuvering.  Burcher 

[8]  discusses  both  model  testing  and  mathematical  approaches  for  predicting  ship 
maneuvering  forces.  Crane  et  al.  [9]  give  another  general  overview  of  ship  maneuvering 
forces.  Inoue  et  al.  [10]  present  equations  for  estimating  hull  maneuvering  forces  based 
on  model  test  data. 

When  considering  the  maneuvering  literature  in  the  context  of  ShipMo3D,  it  should 
be  noted  that  ShipMo3D  has  been  developed  based  primarily  on  the  theory  for  ship 
seakeeping,  which  uses  somewhat  different  definitions  and  approaches  from  those  in 
the  ship  maneuvering  held.  The  first  major  difference  between  the  maneuvering  and 
seakeeping  approaches  is  that  maneuvering  predictions  are  usually  based  on  ship-fixed 
axes,  while  seakeeping  predictions  are  usually  based  on  translating  earth  axes.  The 
second  major  difference  is  that  maneuvering  predictions  usually  assume  that  hull 
forces  are  independent  of  frequency  of  motion,  while  seakeeping  predictions  duely 
consider  the  influence  of  motion  frequency.  Consequently,  care  must  be  taken  when 
applying  results  from  the  maneuvering  literature  to  ShipMo3D. 

3  Axis  Systems  and  Equations  of  Motion  for 
ShipMo3D  Computations 


Previous  ShipMo3D  reports  [1,  2,  3,  4]  consider  the  motions  in  waves  of  a  ship  with 
quasi-steady  speed  and  heading.  This  report  introduces  motions  for  a  freely  maneu¬ 
vering  ship,  which  must  ultimately  be  expressed  using  earth-fixed  axes.  ShipMo3D’s 
earth-fixed  axis  system  is  shown  in  Figure  1,  and  was  first  introduced  in  Reference  3. 
The  instantaneous  heading  (to)  of  the  ship  is  denoted  x  and  the  heading  (from)  of 
ocean  waves  is  denoted  v. 

In  keeping  with  ShipMo3D’s  seakeeping  pedigree,  the  ship  motions  are  solved  using 
the  translating  earth  coordinate  system  shown  in  Figure  2.  When  ocean  waves  are 
present,  their  direction  relative  to  the  ship  is  given  according  to  Figure  3.  For  large 
amplitude  motions,  forces  on  the  ship  can  be  evaluated  using  the  ship  based  axis 
system  of  Figure  4  and  then  transformed  to  translating  earth  axes. 

When  computing  ship  motions  in  the  time  domain,  it  is  necessary  to  select  a  suitable 
nominal  heading  and  nominal  speed  for  orienting  the  translating  earth  coordinate  sys- 
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tem  used  to  solve  the  ship  motions  at  each  time  step.  Computations  with  ShipMo3D 
indicate  that  the  translating  earth  coordinate  system  can  be  oriented  based  on  the 
instantaneous  heading  and  ship  speed  along  that  heading  at  each  time  step. 

During  each  time  step  of  a  time  domain  simulation,  the  forces  and  resulting  ac¬ 
celerations  are  computed  in  translating  earth  axes.  The  resulting  accelerations  in 
earth-fixed  axes  are: 


Xf  =  f/i  cos  x  ~  >h  sin  x 
VS  =  ~Vi  sin  x  +  V2  X 

x  = 


(1) 

(2) 
(3) 


Heave,  roll,  and  pitch  are  equivalent  in  earth-fixed  and  translating  axis  systems. 

The  equations  of  motion  in  the  time  domain  include  terms  presented  in  earlier  re¬ 
ports  [2,  4],  and  new  terms  introduced  in  this  report  that  are  important  for  a  freely 
maneuvering  ship.  These  new  terms  include  forces  due  to  resistance,  propulsion, 
rudder-propeller  interaction,  hull  maneuvering  (which  can  include  lift  and  viscous 
contributions),  and  hull  cross- flow  drag.  The  equations  of  motion  are: 

([M]  +  [H(C,oo)]) 

-  B(U,oo,  if) 

-  C(U,  oo,  rj ) 

+  {Fresist(t)} 

rudder— prop 

The  ship  inertia  matrix  [M]  is  based  on  the  dry  ship.  The  added  mass  matrix 
[A{U,  ex:)]  is  the  infinite  frequency  value  dependent  on  ship  speed  U,  and  includes 
contributions  from  the  hull  and  appendages  as  follows: 

Nfoil  N rudder 

[A(U,  oo)]  =  [Ahull(U,  oo)]  +  X  [■ 4 +  XI  [Arudder~i] 

i—  1  i=  1 

Nskeg  Nm 

+  X  [A3^9^]  +  X  [ A bk~{] 

i=l  i= 1 

where  {Ahull\  is  the  hull  added  mass,  is  the  added  mass  for  static  foil  i  (e.g., 

propeller  shaft  bracket),  [Arudder~l]  is  the  added  mass  for  rudder  i,  [. Aske9~l ]  is  the 
added  mass  for  skeg  i,  and  [Abk~l]  is  the  added  mass  for  bilge  keel  i.  The  appendage 
added  masses  are  independent  of  speed  and  frequency  of  motion.  The  damping  matrix 


(5) 


W)}  = 

W)}  -  S'  [V',”“(C/,«-t)]{'I)(t)}  dr 

J  — OO 

{<?«}  +  +  {FD(t)}  +  {F-»"(()} 

_|_  |_pPr°P(t)|  -)-  {^rudder-deflect 

(t)} 


(4) 
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includes  contributions  from  the  hull  and  appendages  as  follows: 


B(U,  00,77) 


j^hull—mannj  ^ 

N foil 


^phull— rad^jj^  j  _|_ 

+  [Bhull~visc(U,  |t)4|)]  +  ^  [Bfoil-\U,  1 774 1 ) ] 

2=1 

^^r  udder  skeg 

+  [Brudder~,(u,  |ij4|)]  +  £  [B“'9-‘(t/,|,i4|)] 

2=1 
Nbk 

+  £  [B“"(C  |^D] 


2=1 


(6) 


2=1 


where  is  the  hull  frequency  independent  damping  clue  to  the  potential  flow, 

^hull-man j  jg  j-jie  }m]|  maneuvei'ing  damping,  which  can  include  forces  from  lift  and 
cross-flow  drag,  [Bhull~msc ]  is  the  hull  viscous  damping  due  to  roll,  ^Byml~t]  is  the 
damping  for  static  foil  i,  [. Brudder~l ]  is  the  damping  for  rudder  i.  is  the 

damping  for  skeg  i,  and  [Bbk~l]  is  the  damping  for  bilge  keel  i.  The  retardation 
function  matrix  [Khul1]  is  due  to  the  dependence  of  hull  hydrodynamic  forces  on 
frequency  of  motion.  The  stiffness  matrix  includes  contributions  from  buoyancy,  hull 
potential  flow  forces,  and  hull  and  appendage  lift  forces  as  follows: 


C(U,  00, 77)  =  [Cbuoy]  +  [chull{U)\  + 


C 


hull— man 


(Ujl) 


N 


foil 


N rudder 


Y  [c^-fu^  +  Y  [c^^-fu)} 

2=1  2=1 

N skeg  Nbk 

Y  [Cske9~l(u)\  +  J]  [Cbk-\U)] 


(7) 


2=1 


2=1 


where  [Cbuoy]  is  the  hull  buoyancy  stiffness,  [chul1  is  the  hull  frequency  independent 
stiffness  clue  to  the  potential  flow,  ^ChuU~man]  is  the  hull  stiffness  due  to  maneuvering 
forces,  [CFd~l]  is  the  lift  stiffness  for  static  foil  i,  [Crudder~l]  is  the  lift  stiffness  for 
rudder  i,  [Csfce9_*]  is  the  lift  stiffness  for  skeg  i,  and  [Cbk~l]  is  the  lift  stiffness  for  bilge 
keel  i.  The  incident  wave  excitation  vector  {F1}  consists  of  the  following: 

^ f  oil  Nruc[cier 

{F\t)}  =  {FI~hull(t)}  +  Y  +  Y  {FI-rudder~i{t)} 

2=1  2=1 
Nskeg  Nbk 

+  {F^^^ft)}  +  (8) 

2=1  2=1 

where  {FI-hullj  is  the  incident  wave  force  on  the  hull,  is  the  incident 

wave  force  on  foil  i,  ^BI~rudder~'1}  is  the  incident  wave  force  on  rudder  i,  |F/-sfce9-*} 
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is  the  incident  wave  force  on  skeg  i,  and  ^F1^bk^'1}  is  the  incident  wave  force  on  bilge 
keel  i.  The  wave  diffraction  force  vector  {FD}  typically  only  consists  of  hull  forces, 
unless  panelled  appendages  are  included  in  diffraction  computations  [4], 

The  hull  cross-flow  drag  { Fcross } ,  ship  resistance  { Fresist },  propulsion  {. Fprop },  rudder 
deflection  ^Frudder~de^lect} }  and  rudder-propeller  interaction  {. Frudder~prop }  are  impor¬ 
tant  force  components  for  a  freely  maneuvering  ship,  and  are  discussed  in  Sections  5 
to  9. 

Some  ship  force  components  are  most  easily  computed  using  ship-fixed  axes,  followed 
by  translation  to  translating  earth  axes.  Figure  4  shows  ship-based  axes,  which 
can  vary  significantly  from  translating  earth  axes  when  motion  amplitudes  are  large. 
Within  ShipMo3D,  ship  resistance  and  propulsion  forces  are  assumed  to  act  along  the 
xs  direction  of  ship-based  axes.  If  buoyancy  and  incident  wave  forces  are  computed 
using  integration  of  panel  pressures  in  the  time  domain,  then  these  components  are 
also  evaluated  using  ship-fixed  axes. 

After  selected  force  components  are  computed  in  ship-based  axes,  they  are  trans¬ 
formed  to  translating  earth  axes  using  the  following: 


W) 

=  F,  (f)  COS  7/5  COS  776 

—  F2  (t)  (—cos r/4  sin  r/6  +  sin 7/4  sin 7/5  cos tjq) 

+  F3  (t)  (cos  r/4  sin  7/5  cos '7/6  +  sin  7/4  sin  t/6) 

(9) 

W) 

=  Ff(t)  COS  '7/5  sin  7/6 

+  Ff  (t)  (cos  7/4  cos  7/6  +  sin  7/4  sin  7/5  sin7/6) 

+  F3  (t)  (cos  7/4  sin  7/5  sin  7/6  —  sin  7/4  cos  r/6) 

(10) 

W) 

=  —  Ff(t)  sin  7/5  +  F29(f)  sin  7/4  cos  7/5 

+  F3  (t)  cos  7/4  cos  7/5 

(11) 

W) 

=  Ff(t) 

(12) 

W) 

=  F55  (t)  cos  7/4  -  F^f  (t)  sin  7/4 

(13) 

W) 

=  F55  (t)  sin  7/4  +  F65  (f )  cos  7/4 

(14) 

where  Fj  is  the  force  for  mode  j  in  translating  earth  axes,  and  F-  is  the  force  for 
mode  j  in  ship-based  earth  axes. 

For  a  ship  executing  a  turning  circle,  several  terms  are  used  to  describe  the  turning 
circle  characteristics,  as  illustrated  in  Figure  5.  Dimensions  are  typically  referenced 
to  the  point  at  which  rudder  deflection  was  initiated.  The  advance  refers  to  the 
displacement  in  the  direction  of  the  initial  heading,  and  the  transfer  refers  to  the 
displacement  lateral  to  the  the  direction  of  the  initial  heading.  The  terms  advance 
and  transfer  are  both  often  applied  to  either  any  point  on  a  turning  circle  or  to  the 
specific  point  at  which  the  ship  heading  has  changed  by  90  degrees  relative  to  the 
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original  heading.  The  tactical  diameter  is  the  transfer  when  the  heading  has  changed 
by  180  degrees. 


4  Hull  Maneuvering  Forces 


To  provide  accurate  maneuvering  predictions,  ShipMo3D  requires  new  hull  maneu¬ 
vering  force  terms  that  have  not  been  previously  included  for  ShipMo3D  seakeeping 
calculations.  These  new  hull  force  terms  are  due  primarily  to  forces  arising  from 
lift  on  the  hull  and  cross-flow  drag.  The  hull  maneuvering  force  terms  are  expressed 
within  ShipMo3D  such  that  they  can  be  determined  from  standard  maneuvering  test 
results  if  available. 


The  prediction  of  hull  maneuvering  forces  is  discussed  in  many  references.  Ship- 
Mo3D  estimates  hull  maneuvering  forces  based  on  Inoue  et  al.  [10],  who  developed 
equations  based  on  physical  considerations  and  analysis  of  model  test  data.  Sway  and 
yaw  damping  coefficients  include  both  linear  and  nonlinear  terms,  and  are  expressed 
as: 


)hull—man 

22 

tt  td  hull— man 

~  U  t>22U 

+ 

u 

\v' 

|  ryhull—man 

1  n22U\v'\ 

+ 

U  |  r' 

|  jDhull—man 

1  n22U\r'\ 

(15) 

>hull—man 

26 

tt  td hull— man 

~  U  ^26U 

+ 

u 

\r' 

|  jDhull—man 

1  n26U\r'\ 

(16) 

ihull—man 

62 

tt  m hull— man 

—  U  -^62 U 

+ 

u 

r'2 

TDhull—man 

n62Ur'2 

(17) 

► hull— man 
66 

tt  td  hull— man 

+ 

u 

\r' 

|  -jDhull—man 

1  n66U\r'\ 

+ 

U  v’2 

jDhull—man 

n66Uv'2 

(18) 

The  non-dimensional  sway  velocity  v'  is  given  by: 


The  non-dimensional  yaw  rate  is  given  by: 

/  he  L 

r  =  - 

U 


(19) 


(20) 


Inoue  et  al.  present  methods  for  computing  the  relevant  terms  in  Equations  (15)  to 
(18).  When  computing  linear  terms  for  the  above  equations,  it  must  be  noted  that 
equations  in  Inoue  et  al.  include  contributions  due  to  hull  radiation  forces,  which 
are  already  included  in  the  ShipMo3D  library.  The  resulting  hull  maneuvering  terms 
computed  using  nondimensional  coefficients  based  on  Inoue  et  al.  are: 

\  P  L  TmU  y;  -  =  oo)  (21) 

\p  £2  Tmid  y;  -  BS~™d(u,e  =  oo)  (22) 


jDhull—man 

n22U 


td  hull— man 
^26 
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Advance  (|Ax|  =  90°) 


-  Tactical  diameter  (|Ax|  =  180°) 

Transfer  (|Ax|  =  90°)  — >\ 


Figure  5:  Turning  Circle  Characteristics 
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BhuU-man  =  Tmid  =  Oo)  (23) 

BhuU-man  =  ApL3Tm.dN'_BhnU-rad{LUe  =  oo)  (24) 


where  Tmid  is  draft  at  midships.  The  hull  radiation  damping  terms  to  be  subtracted 
from  the  above  equations  are  given  in  Reference  2  as  follows: 


B ? 


hull— rad . 


j2U 


UP  =  oo  = 


r>hull—rad  /  ,  \ 

Bj6U  K  =  OO)  = 


' — "  hull  , 

<9^2  (t/  =  0,u;e  =  oo) 
dx 

-A)f\U  =  0,cue  =  oo) 

' — "  /m/Z 

(9A,6  (17  =  0,o;e  =  oo) 


(9a: 


(25) 


(26) 


where  <9/1  /dx  is  determined  from  the  x  derivative  of  potentials  giving  added  mass 
and  Aj%11  is  added  mass  due  to  sway  motion,  with  U  =  0  denoting  zero  forward  speed 
and  u)e  —  oo  denoting  infinite  encounter  frequency. 


The  linear  nondimensional  coefficients  for  hull  maneuvering  forces  are  obtained  based 
on  Inoue  et  ah,  with  sign  changes  to  conform  to  conventions  commonly  used  within 
the  maneuvering  literature: 


y:  = 


y:  =  - 


N' 

r 


T, 


IT 


mid 


+  1-4  Cb  -jt 


N’  = 


Tmid 

2  T~ 
2  Tmid 
L 

0.54 


1  +  0.8 


t 


stern 


1  - 

2  Tmid 


2T, 


mid 


l  + 


T7 

mid  , 

0-27  tatern\ 

la  Tmid  ) 


2  t stern  \ 

3  Tmid  ) 


1  +  0.30 


Astern 


T 


mid 


(27) 

(28) 

(29) 

(30) 


where  Cb  is  block  coefficient,  B  is  beam,  and  tstern  is  trim  by  stern.  The  term  in 
Equation  (29)  is  given  by: 


2  Tmid 

7T  Tmid  +  1-4  Cb  B 


(31) 


The  nonlinear  hull  maneuvering  force  terms 


ryhull—man 

n22U\v'\ 


mhull—man 

n22U\r'\ 


1 

2 
1 

2 


are  given  by: 

P  T  Tmid  Yv\v\ 
P  L  Tmid  Ty|y| 


(32) 

(33) 
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B 

B 

B 

B 


hull— man 
26U\r'\ 

hull— man 
62  Ur'2 

hull— man 
66U\r'\ 

hull— man 
66Uv'2 


~  2  P  L2  Tmid  ^r\r\ 

—  g  P  L2  Tmid  N'vr 2 

—  2  P  Tmid  ^r\r\ 

~  2  P  Tmid  N'rv2 


(34) 

(35) 

(36) 

(37) 


The  non-dimensional  terms  in  the  above  equations  are  evaluated  based  on  Inoue  et 
ah,  with  signs  changed  to  conform  to  conventions  commonly  used  in  the  maneuvering 
literature: 


Y' 

*  v\v\ 


N'  2 

vrz 


NL 


|r| 


N' 


0.09  6.5  (1  CB)  TmJf 

JD 

(38) 

0.44  +  1.78  (1  CB)  TmJf 

JD 

(39) 

0.0 

(40) 

0.0 

(41) 

(  -0.060 

for 

clB  <  0.06 

<  0.146 +  1.8C®S  6(cfB)2 

for 

0.06  <  Cblb  <  0.2 

(42) 

[  -0.026 

for 

CbB  >  0.2 

-0.2 

(43) 

The  terms  YY  and  NY2  in  the  above  equations  are  set  to  zero  because  of  the  signif¬ 
icant  scatter  in  measured  values. 


Experimental  data  presented  by  Inoue  et  ah  suggest  that  actual  nondimensional  ma¬ 
neuvering  coefficients  can  vary  significantly  from  those  predicted  using  Equations 
(27)  to  (30)  and  (38)  to  (43).  Table  1  gives  nominal  uncertainties  associated  with 
estimated  hull  force  coefficients.  The  uncertainties  in  Table  1  have  been  set  such  that 
actual  maneuvering  coefficients  arc  likely  to  fall  within  the  specified  range. 


The  maneuvering  forces  presented  thus  far  consider  only  sway  and  yaw.  In  the  current 
work,  the  vertical  location  of  the  lateral  hull  maneuvering  force  is  assumed  to  act  at 
the  waterline  based  on  Bertram  [5].  The  location  of  the  hull  maneuvering  force 
relative  to  the  centre  of  gravity  is  given  by: 

—hull— man  =  _Z<JG  (44) 

where  z^f  is  the  location  of  the  ship  centre  of  gravity  relative  to  the  calm  waterline. 
The  above  vertical  location  is  also  assumed  to  be  the  reference  point  for  lateral 
motions  influencing  hull  maneuvering  forces.  The  following  terms  can  then  be  derived: 


m hull— man  _  _  td hull— man  —hull— man 

_d24  —  D  22  Z 


(45) 
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T)  hull— man 
-^42 

TD  hull— man 
~  ^22 

—hull— man 

(46) 

td  hull— man 
-O44 

_  hull— man  ^ 

—hull—man  \  ^ 

(47) 

td  hull— man 
-^46 

TD  hull— man 
~  ^26 

—hull— man 

(48) 

TD  hull— man 
-^64 

TD  hull— man 
~  -^62 

—hull— man 

(49) 

It  should  be  noted  that  the  approach  developed  by  Inoue  et  al.  for  predicting  hull 
maneuvering  coefficients  is  intended  for  hulls  with  lateral  profiles  approximating  a 
rectangle.  For  ships  such  as  frigates  having  significant  variation  of  draft  with  lon¬ 
gitudinal  location,  care  should  be  taken  in  applying  the  method  of  Inoue  et  al.  for 
predicting  hydrodynamic  coefficients. 

Table  1:  Uncertainties  for  Estimated  Nondimensional  Hull  Maneuvering  Coefficients 


Y' 

V 

±0.05 

Y' 

r 

±0.05 

N' 

1  V 

±0.02 

N'r 

±0.02 

Y' 

*v\v\ 

±0.2 

Y' 

*  v\r\ 

±0.2 

nVi 

±0.05 

Kr> 

±0.05 

Nr\r\ 

±0.02 

N'rv* 

±0.2 

5  Separate  Treatment  of  Hull  Cross-Flow 
Drag  Forces 


The  hull  maneuvering  force  terms  presented  in  the  previous  section  typically  include 
the  effects  of  cross-flow  drag.  Alternatively,  it  is  possible  evaluate  hull  cross-flow  drag 
separately,  as  discussed  by  Bertram  [5].  For  a  ship  exposed  to  a  uniform  cross-flow 
velocity  of  vcross ,  the  cross-flow  drag  force  will  be: 

Fcross  =  ^  pvcross  \vcros8\  Ay  CDy  (50) 

where  Ay  is  the  hull  lateral  profile  area.  The  hull  cross-flow  drag  velocity  acting  on  a 
maneuvering  ship  is  caused  by  the  sway  and  yaw  motion  of  the  ship,  and  varies  along 
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the  ship  hull.  If  cross-flow  drag  is  computed  separately  from  other  hull  maneuver¬ 
ing  forces  within  ShipMo3D,  then  the  hull  is  divided  into  a  number  of  longitudinal 
segments  (typically  10-20)  for  evaluation  of  the  cross-flow  drag  as  follows: 

NSeg 

Fr„  =  _  p  E  \vcro,,(xAs-i)\  V.  (51) 

1 

where  Nseg  is  the  number  of  hull  longitudinal  sections,  XAy-i  is  the  nominal  x  coordi¬ 
nate  of  section  i,  Ay_i  is  the  lateral  profile  area  of  section  i,  and  Cr>y  is  the  cross-flow 
drag  coefficient.  The  local  hull  cross-flow  velocity  is  given  by: 


v"°”(XAy-i)  =  -1)2  -  XAV-i  1(6 


(52) 


The  roll  and  yaw  cross-flow  drag  forces  are  evaluated  in  a  manner  similar  to  the  sway 
cross-flow  drag  force: 

Nseg 

FZ°SS  =  \vCrOSS(XAy-i)\  Ay-iZAy-iCDy  (53) 

i= 1 
Nseg 

FZ°SS  =  -P  ^^{XAy-i)  \vCr°88(XAy-i)\  Ay-i  XAy-i  CDy  (54) 

i= 1 

where  ZAy-i  is  the  nominal  z  coordinate  of  the  underwater  prohle  for  hull  longitudinal 
segment  i.  The  geometric  properties  of  the  hull  longitudinal  segments  are  evaluated 
based  on  the  wet  panelled  hull  surface  using  the  following: 


A  — 

F^y-i 

Akij'ij )  |  fly— k(i,j)  | 

3= 1 

Nv-i 

%Ay—i 

1 

JJ  L 

^  ^  Akij'i  j )  \Hy— k(i,j) 

3= 1 

|  Xk(i,j ) 

9  A 

^^i.y  —  1 

Nv-i 

%Ay—i 

1 

JJ  L 

£  Ak{i,j)  |ny-fc(*,i)| 

3= 1 

Zk(i,j) 

9  A 

(55) 

(56) 

(57) 


where  Np_l  is  the  number  of  hull  panels  whose  centroid  lies  within  hull  longitudinal 
segment  i,  k(i,j )  is  the  global  hull  panel  index  for  panel  j  on  longitudinal  segment  i, 
Ak  is  the  area  of  hull  panel  k ,  ny_k  is  the  y  normal  component  for  hull  panel  k ,  x k 
is  the  x  coordinate  of  the  centroid  of  hull  panel  k,  and  Zk  is  the  z  coordinate  of  the 
centroid  of  hull  panel  k. 


Bertram  presents  lateral  drag  coefficients  ranging  from  0.56  to  0.98  for  four  different 
ships.  It  is  expected  that  lateral  drag  coefficients  for  most  ships  will  be  between  0.5 
and  1.1. 
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As  previously  noted,  hull  cross-flow  drag  forces  are  typically  included  in  hull  maneu¬ 
vering  force  coefficients.  When  including  cross-flow  drag  in  maneuvering  coefficients, 
there  appears  to  be  an  implicit  assumption  that  the  cross-flow  velocity  is  small  rela¬ 
tive  to  ship  forward  velocity.  In  some  cases,  such  as  towing  operations,  the  cross-flow 
velocity  can  be  a  significant  fraction  of  the  ship  forward  velocity.  In  such  cases,  it 
is  likely  more  appropriate  to  model  cross-flow  drag  forces  separately  from  other  hull 
maneuvering  forces. 


6  Hull  Resistance  Forces 


Several  references  provide  good  overviews  on  hull  resistance  forces,  including  Bertram 
[5],  Newman  [11],  and  van  Manen  and  van  Oossanen  [12].  ShipMo3D  evaluates  hull 
resistance  forces  using  input  resistance  coefficients  which  are  given  as  a  function  of 
ship  speed: 

F[esist  =  — l-pU\U\AwCDx{U )  (58) 

where  Aw  is  the  hull  wetted  surface  area  and  Cdx  is  the  ship  resistance  coefficient. 
Experimental  data  and  CFD  computations  are  among  the  available  methods  for  se¬ 
lecting  suitable  resistance  coefficient  values. 

7  Propeller  Forces 


Modelling  of  propeller  forces  is  covered  in  many  references,  including  Holtrop  and 
Mennen  [13,  14],  Holtrop  [15],  and  Van  Manen  and  Van  Oossanen  [16].  The  thrust 
created  by  the  propeller  can  be  expressed  by: 


FPr°P  =  (1  -  W)  P  nprop  ytop  KT(Jprop)  (59) 


where  tprop  is  the  propeller  thrust  deduction  coefficient,  nprop  is  the  propeller  speed 
in  revolutions  per  second,  and  Dprop  is  the  propeller  diameter.  The  propeller  thrust 
coefficient  Kt  is  a  function  of  the  advance  coefficient  Jprop)  which  is  given  by: 


J, 


prop 


U  (1  Wprop ) 

^ prop  Dpr0p 


(60) 


where  wprop  is  the  propeller  wake  fraction.  Both  the  thrust  deduction  coefficient  fprop 
and  propeller  wake  fraction  wprop  tend  to  increase  with  hull  block  coefficient,  and  are 
generally  larger  for  single  screw  ships  than  for  twin-screw  ships.  Thrust  deduction 
coefficients  typically  fall  within  the  range  0.0  <  tprop  <  0.2,  and  propeller  wake 
fractions  typically  fall  within  the  range  0.0  <  wprop  <  0.4. 
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Figure  6:  Example  Curve  of  Propeller  Thrust  Coefficient  Kt  Versus  Advance  Coef¬ 
ficient  Jprop  (adapted  from  Reference  16) 

The  propeller  thrust  coefficient  Kt  decreases  with  increasing  advance  coefficient  Jprop- 
Figure  6  shows  an  example  thrust  coefficient  curve  adapted  from  Reference  16. 

ShipMo3D  assumes  that  the  propeller  thrust  acts  along  the  ship  longitudinal  axis. 
Transformation  to  force  components  in  translating  earth  axes  is  done  using  Equations 
(9),  (10)  and  (11). 
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8  Rudder  Deflection  Forces 


Forces  arising  from  deflection  of  rudders  (and  other  deflectable  appendages)  are  com¬ 
puted  as  follows: 


{F 


rudder— deflect 


(«)}  =  -  E  [cy^do]  S 

i= 1 

N rudder 

|  ^  ^  prudder— drag— i  ^ 

2—1 


rudder— i 


(61) 


where  [C'^cWer  *]  is  the  stiffness  matrix  for  ship  forces  arising  from  rudder  deflections 
for  rudder  i,  $rudder~l  is  the  deflection  of  rudder  i,  and  ^Frudder~dra9^1}  js  the  drag 
forces  vector  on  rudder  i.  ShipMo3D  uses  a  rudder  deflection  convention  of  positive  for 
counter-clockwise  rotation  when  viewed  from  inside  the  ship.  The  rudder  deflection 
stiffness  terms  are  evaluated  based  on  Reference  4.  ShipMo3D  calculations  for  rudder 
drag  forces  are  based  on  fnoue  et  al.  [17]  as  follows: 

pr udder— deflect  prudder  ^ru^der  (62) 

where  F™‘dd' '  is  the  total  force  normal  to  the  rudder,  including  contributions  from 
ship  motions  in  waves.  The  force  normal  to  the  rudder  is  evaluated  as  follows: 


r>  rudder 


F. 


rudder 


rprudder 

N 


6 

6 

6 

\  A  prudder 
/  -j 

Vj 

\  A  sv rudder 

-  Vj 

\  '  s~i rudder  crudder 

~  t'25  ° 

(63) 

3= 1 

3= 1 

3=1 

6 

6 

6 

\  A  prudder 
/  - j  3j 

Vj 

\  ^  sy rudder 

~  °3 j  Vj 

^  ^  pirudder  ^ rudder 

(64) 

3= 1 

3  =  1 

3= 1 

pruMe.  cos  Trudder  -  F™dder  sin  r rudder  (65) 


where  F^udder  i§  the  total  rudder  lateral  force,  B™dder  is  rudder  force  damping  due  to 
ship  motion  mode  j,  C™dder  is  rudder  force  stiffness  due  to  ship  motion  mode  j,  C^dder 
is  rudder  force  stiffness  due  to  rudder  deflection,  F^udder  is  total  rudder  vertical  force, 
and  Vrudder  is  rudder  dihedral  angle  (0°for  rudder  to  port,  90°for  rudder  upward). 


9  Rudder-Propeller  Interaction  Forces 


The  rudder  forces  presented  in  Reference  4  and  in  the  previous  section  do  not  con¬ 
sider  the  influence  of  propellers.  To  increase  the  steering  efficiency  of  rudders,  they 
are  typically  placed  within  the  slipstream  of  ship  propellers.  The  resulting  rudder 
forces  are  significantly  greater  than  if  the  rudder  were  placed  outside  of  the  propeller 
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slipstream.  ShipMo3D  models  the  influence  of  the  propeller  on  rudder  forces  using 
an  approach  developed  by  Soding  [18].  In  addition  to  the  rudder  lift  proportional  to 
U2  presented  in  Reference  4,  the  presence  of  the  propeller  induces  the  following  lift 
and  drag  on  a  rudder  in  the  propeller  slipstream: 


rpr  udder— prop 
"lift 

j-ir  udder— prop 
drag 


pirudder—prop  pprop  I  \  - |_ 


pi  rudder— prop  pprop  /  ^ 


\/l  +  Cth 
1 


th 


Sin  judder 

)  (1  -  cos  Srudder) 


(66) 

(67) 


where  Crudder~prop  is  a  rudder-propeller  interaction  coefficient,  and  Ct,h  is  the  thrust 
loading  coefficient  of  the  propeller.  The  rudder- propeller  interaction  coefficient  nor¬ 
mally  ranges  between  0.0  and  1.0,  increasing  with  proximity  of  the  rudder  to  the 
propeller  and  with  size  of  the  rudder  relative  to  the  propeller.  The  effect  of  increas¬ 
ing  rudder  size  becomes  negligible  when  the  rudder  span  exceeds  twice  the  propeller 
diameter.  The  propeller  thrust  loading  coefficient  Cth  is  given  by: 


Cth  = 


pprop 


1/2  p  f/2  (1  -  wpropy  vr/4  Dl 


prop 


(68) 


The  rudder-propeller  interaction  forces  in  ship-based  axes  are  as  follows: 


jp  rudder— prop,  S  _ 

rprudder—prop 

drag 

(69) 

T^r  udder— prop,  S 

T~\r udder— prop  •  -n 

~Fuft  smr 

(70) 

T^r  udder— prop,  S 

rprudder—prop  -n 

Fuft  cos  r 

(71) 

rpr  udder— prop,  S 

"a  — 

— S  jpr udder— prop, S  —S  jpr udder— prop,  S 

y  3  %  2 

(72) 

T^r  udder— prop,  S 
^5  — 

-  (xs  +  c/4)  Fludder-prop's 

|  —S  pr udder— prop,  S 

(73) 

rpr  udder— prop,  S 

[xS  +  c/4)  F™dder-prop,S  _ 

— S  pr  udder— prop,  S 

y  i 

(74) 

where  xs,ys,zs  is  the  location  of  the  rudder  centroid  in  ship-based  axes,  and  c  is 
the  mean  rudder  chord  length.  Conversion  from  ship-based  axes  to  translating-earth 
axes  is  done  using  Equations  (9)  to  (14). 


When  modelling  rudder-propeller  interaction  forces,  the  greatest  uncertainty  lies  in 
the  selection  of  the  rudder-propeller  interaction  coefficient  Crudder~prop .  If  a  large 
rudder  is  located  immediately  aft  of  a  propeller,  then  the  interaction  coefficient 
Cr udder -prop  can  set  i  q  Alternatively,  if  the  rudder  is  clearly  not  in  the  pro¬ 
peller  slipstream,  then  the  interaction  coefficient  crudder~prop  can  be  set  to  0.0.  For 
intermediate  cases  such  as  typical  modern  frigates  with  a  single  rudder  aft  of  two 
propellers,  it  is  more  difficult  to  select  a  suitable  interaction  coefficient. 
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10  Rudder  and  Propeller  RPM  Control 
Modelling 


Simulation  of  the  motions  of  a  freely  maneuvering  ship  requires  modelling  of  rudder 
and  propeller  RPM  controls.  ShipMo3D  uses  linear  control  systems  to  model  both 
rudder  and  propeller  responses  to  input  commands. 

The  rudder  control  system  is  the  same  as  that  presented  in  Reference  4  for  a  ship 
travelling  with  quasi-steady  speed  and  heading.  The  rudder  response  characteristics 
are  modelled  as  follows: 


£ rudder 


.  o  >  rudder 

+  ^  Qs  VS 


g rudder 


i  ,  ,2  rr udder 
+  UJS  d 


,  ,2  rr udder 

US  (>c 


(75) 


where  Srudder  is  rudder  acceleration,  (s  is  the  nondimensional  damping  response  con¬ 
stant,  is  the  rudder  response  natural  frequency,  Srudder  is  rudder  velocity,  and 
^rudder  commanc[  rudder  angle.  The  rudder  motions  in  the  time  domain  are  easily 
evaluated  as  follows: 


cr udder  ,  ,2  (  rr udder  cr udder \  r>  /•  ,  ,  rr udder 

d  =  us  (dc  -d  J  -  2  C sVfid 


(76) 


The  autopilot  for  a  freely  maneuvering  ship  is  very  similar  to  the  autopilot  for  a  ship 
with  quasi-steady  speed  and  heading  presented  in  Reference  4.  The  main  difference 
is  that  an  autopilot  for  a  freely  maneuvering  ship  uses  input  motions  in  earth-fixed 
axes,  while  an  autopilot  for  a  ship  with  quasi-steady  speed  and  heading  uses  input 
motions  in  translating  earth  axes.  Consequently,  yaw  gains  for  a  freely  maneuvering 
autopilot  are  of  opposite  sign  to  those  for  an  autopilot  based  on  translating  earth 
axes. 


The  command  rudder  angle  for  a  freely  maneuvering  autopilot  is  given  by: 

=  E  [n,  a’  +  vSi  M  (,,/  -  4)' 

3= 1 


(77) 


where  is  the  acceleration  gain  for  motion  mode  j,  is  the  velocity  gain,  kdj  is 
the  displacement  gain,  and  rfC]  is  the  command  displacement. 

The  simulation  of  a  rudder  control  system  including  autopilot  requires  selection  of 
suitable  input  values.  For  a  US  Coast  Guard  cutter  representative  of  modern  frigate 
design,  Smith  [19]  indicates  a  maximum  rudder  deflection  of  35  degrees  and  maximum 
rudder  rate  of  3  degrees  per  second.  For  modelling  of  a  conventional  downward 
rudder  using  ShipMo3D,  the  yaw  gain  and  yaw  velocity  gain  for  a  freely  maneuvering 
autopilot  will  typically  have  values  less  than  zero. 


18 


DRDC  Atlantic  TM  2005-071 


The  propeller  control  system  is  modelled  very  similarly  to  the  rudder  control  system, 
with  the  rate  of  change  of  propeller  RPM  being  modelled  as  follows: 


RPMprop  =  u2rpm  {RPMq°p  —  RPMprop)  —  2  Crpm  ^rpm  RPMpr°P  (78) 

where  RPM  is  the  second  derivative  with  respect  to  time  of  propeller  RPM, 
ujrpm  is  the  RPM  response  natural  frequency,  RPM^fop  is  command  propeller  RPM, 
RPMprop  is  propeller  RPM,  Crpm  is  RPM  response  damping,  and  RPMpr°P  is  the 
first  derivative  with  respect  to  time  of  propeller  RPM. 

1 1  Validation  of  Turning  Circles  for  Esso 
Osaka  Tanker 


Initial  validation  of  ShipMo3D  maneuvering  predictions  has  been  done  using  turning 
circles  for  the  oil  tanker  Esso  Osaka,  with  particulars  as  given  in  Table  2.  Crane  [20] 
describes  a  comprehensive  set  of  maneuvering  trials  conducted  for  the  Esso  Osaka. 
The  high  quality  of  the  reported  trials  data  has  prompted  much  subsequent  work  on 
prediction  of  maneuvering  forces  and  motions  for  the  Esso  Osaka.  A  recent  Special 
Committee  of  the  International  Towing  Tank  Conference  [6]  has  produced  an  excellent 
overview  on  maneuvering  forces  and  motions  for  the  Esso  Osaka. 

Table  2:  Main  Particulars  for  the  Esso  Osaka  During  Maneuvering  Trials 

Length,  L  (between  perpendiculars)  325  m 

Beam,  B  53  m 

Midships  draft,  Tmui  21.73  m 

Trim  by  stern,  tstern  0.0  m 

Displacement,  A  319,400  tonnes 

Longitudinal  CG  forward  of  midship  10.3  m 


Esso  Osaka  geometric  data  provided  by  Crane  were  used  to  develop  ShipMo3D  input. 
Detailed  resistance  and  propulsion  data  were  not  found;  thus,  it  was  necessary  to 
estimate  these  characteristics.  Resistance  coefficients  were  taken  from  the  Series  60 
hull  with  Cr  =  0.80  [21],  which  has  geometric  properties  very  similar  to  the  Esso 
Osaka  (see  Table  3).  Table  4  shows  the  resulting  estimated  resistance  coefficients 
scaled  to  speeds  for  the  Esso  Osaka, 

For  propulsion  characteristics  of  the  Esso  Osaka,  Crane  indicates  that  the  propeller 
diameter  Dvrop  is  9.1  m  and  that  the  propeller  rotates  at  51  RPM  when  the  ship 
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travels  at  10  knots  in  calm  water.  Based  on  data  given  in  Reference  16,  the  wake 
fraction  wprop  was  estimated  to  be  0.352  and  the  thrust  deduction  coefficient  tprop 
was  estimated  to  be  0.2.  To  obtain  a  thrust  coefficient  curve  for  the  Esso  Osaka,  it 
was  decided  that  the  sample  thrust  coefficient  curve  from  Figure  6  could  be  scaled  to 
provide  a  thrust  equivalent  to  the  modelled  resistance  of  the  Esso  Osaka  at  10  knots. 
Surprisingly,  the  sample  thrust  coefficient  curve  provided  a  thrust  within  2  percent 
of  the  required  value,  with  the  final  scaled  thrust  coefficient  curve  for  the  Esso  Osaka 
being  as  follows: 


Kr 


0.394  -  0.197  Jprop 


0.148  Jprop 


(79) 


Table  5  gives  rudder  dimensions  and  modelled  control  properties  for  the  Esso  Os¬ 
aka.  The  rudder  span  and  chord  are  taken  from  Reference  20.  With  the  exception 
of  maximum  rudder  deflection  the  rudder  control  parameters  are  estimates 

based  on  published  values  for  other  ships.  The  rudder-propeller  interaction  coefficient 
(jrudder-prop  estimated  to  be  0.9  due  to  the  proximity  and  large  size  of  the  rudder 
relative  to  the  propeller. 

Figures  7  to  12  present  predicted  and  measured  turning  circle  data  for  the  Esso  Osaka 
with  initial  velocities  of  10  knots  and  7.7  knots.  The  ShipMo3D  predictions  use  hull 
maneuvering  coefficients  based  on  fnoue  et  al.  [10].  The  predictions  show  generally 
good  agreement  with  the  sea  trials  data.  The  most  noticeable  differences  are  that  the 
predicted  ship  speeds  and  absolute  yaw  rates  are  larger  than  actual  values  when  the 
ship  reaches  its  steady-state  turning  condition.  When  examining  the  sea  trial  data, 
the  data  for  an  initial  ship  speed  of  7.7  knots  are  less  consistent  than  the  data  for  10 
knots.  The  lesser  consistency  at  lower  ship  speed  is  likely  due  to  the  ship  being  more 
prone  to  environmental  forces  such  as  those  from  winds  and  currents. 

The  present  ShipMo3D  predictions  for  maneuvering  of  the  Esso  Osaka  are  encour¬ 
aging.  The  final  turning  circle  diameter  for  the  Esso  Osaka  is  equal  to  only  3  ship 
lengths;  thus,  nonlinear  effects  are  significant. 
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Table  3:  Principal  Hull  Coefficients  for  Esso  Osaka  and  Series  60  Hull  with  Cb  =  0.80 


Esso  Osaka 

Series  60  with  C} 

Block  coefficient  Cb 

0.831 

0.80 

Length/beam  L/B 

6.1 

6.5 

Beam/draft  B/T 

2.4 

2.5 

Length  L 

325  m 

122  m 

Table  4:  Estimated  Resistance  Cofffcients  for  Esso  Osaka 


Speed 

(knots) 

2.0 

16.0 

18.0 

19.6 

21.2 

22.9 


Resistance  coefficient 
CDx  =  F/(l/2pU2Aw ) 
0.00306 
0.00306 
0.00308 
0.00328 
0.00370 
0.00426 


Table  5:  Rudder  Dimensions  and  Modelled  Control  Properties  for  Esso  Osaka 


Rudder  span  srudder 

13.85  m 

Rudder  chord  length  crudder 

9.00  m 

Maximum  rudder  deflection  §™dder 

35  degrees 

Maximum  rudder  velocity  4™“er 

3  deg/s 

Response  natural  frequency  u>s 

3  rad/s 

Response  damping  Q 

0.85 

Rudder-propeller  interaction  coefficient  crudder~prop 

0.9 
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Full-scale  trials 
ShipMo3D 


Figure  7:  Turning  Circle  Trajectory  for  Esso  Osaka  Starboard  Turn  at  10  Knots  with 
35  degree  Rudder 


Figure  8:  Yaw  Rate  Versus  Time  for  Esso  Osaka  Starboard  Turn  at  10  Knots  with 
35  degree  Rudder 


Figure  9:  Ship  Speed  Versus  Time  for  Esso  Osaka  Starboard  Turn  at  10  Knots  with 
35  degree  Rudder 
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Figure  10:  Turning  Circle  Trajectory  for  Esso  Osaka  Port  Turn  at  7.7  Knots  with  35 
degree  Rudder 


t  (s) 


Figure  11:  Yaw  Rate  Versus  Time  for  Esso  Osaka  Port  Turn  at  7.7  Knots  with  35 
degree  Rudder 


Figure  12:  Ship  Speed  Versus  Time  for  Esso  Osaka  Port  Turn  at  7.7  Knots  with  35 
degree  Rudder 
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12  Unclassified  Maneuvering  Modelling  of 
HALIFAX 


Examination  of  predicted  maneuvering  forces  and  motions  for  Canada’s  HALIFAX 
class  frigate  can  provide  useful  insight  into  maneuvering  prediction  of  naval  vessels. 
Measured  maneuvering  data  for  HALIFAX  are  classified  and  not  presented  here. 
Similarly,  measured  resistance  and  propulsion  data  for  HALIFAX  are  classified;  thus, 
resistance  and  propulsion  properties  are  modelled  using  approximations  based  on  data 
in  the  open  literature. 

Table  6  gives  main  particulars  for  the  HALIFAX  with  the  loading  condition  based  on 
that  used  for  model  tests  with  the  CPF  hydroelastic  model  [22],  Resistance  coeffi¬ 
cients  given  in  Table  7  have  been  estimated  based  on  values  in  Schmitke  and  Murdey 
[23]  for  a  frigate  with  similar  geometry.  The  tabulated  resistance  coefficients  are 
intended  to  encompass  all  possible  operating  speeds. 

Based  on  the  geometry  of  the  propeller  shaft  brackets  for  HALIFAX,  the  two  HALI¬ 
FAX  propellers  are  modelled  with  propellers  having  nominal  diameters  of  5.6  m.  The 
propeller  wake  fraction  and  thrust  deduction  coefficient  are  both  assumed  to  be  0.0 
based  on  comments  regarding  destroyer  type  vessels  in  Reference  16.  The  propeller 
thrust  coefficient  curve  has  been  designed  based  on  the  following  constraints: 

•  For  a  nominal  ship  speed  of  30  knots,  the  propeller  RPM  should  be  a  suitable 
nominal  value  (e.g.,  100), 

•  The  propeller  thrust  coefficient  curve  is  scaled  from  the  curve  Kt(Jprop)  =  0.40  — 
0.2 JPr0p  —  0.15  Jprop  from  Reference  16, 

•  To  ensure  operation  of  the  propeller  at  a  realistic  advance  ratio,  the  advance 
ratio  at  a  ship  speed  of  35  knots  is  less  than  0.90, 

•  At  zero  advance  coefficient,  the  thrust  coefficient  curve  should  satisfy  0.2  < 

Kt{ Jprop  =  0)  <  0.6. 

The  following  thrust  coefficient  curve  for  the  HALIFAX  propellers  gives  a  nominal 
speed  of  30  knots  for  the  propellers  rotating  at  200  RPM: 

Kt ( Jprop)  =  0.32  —  0.16  Jprop  —  0.12  Jprop  (80) 

Table  8  gives  speeds  in  calm  water  and  associated  propeller  RPM  values  determined 
using  iterative  solution. 

Table  9  gives  modelled  rudder  control  properties  and  the  estimated  rudder-propeller 
interaction  coefficient  for  the  HALIFAX.  With  the  exception  of  maximum  rudder 
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Table  6:  Main  Particulars  for  HALIFAX  Class  Frigate,  CPF  Hydroelastic  Model  Deep 

Departure  Condition 

Length,  L 

124.5  m 

Beam,  B 

14.8  m 

Midships  draft,  Tmid 

4.97  m 

Trim  by  stern,  tstern 

-0.04  m 

Displacement,  A 

4601  tonnes  (salt  water) 

Vertical  centre  of  gravity,  KG 

6.26  m 

Dry  roll  radius  of  gyration  rxx 

5.82  m 

Dry  pitch  radius  of  gyration  ryy 

28.8  m 

Dry  yaw  radius  of  gyration  rzz 

28.8  m 

Table  7:  Resistance  Coefficients  for  Unclassified  Modelling  of  HALIFAX  Maneuvering 

Speed  Resistance  Coefficient 

(knots)  Cr  - 

=  FR/(l/2pU2Aw ) 

<  10 

0.0040 

15 

0.0042 

20 

0.0057 

25 

0.0069 

30 

0.0099 

35 

0.0108 

Table  8:  Ship  Speeds  in  Calm  Water  and  Associated  Propeller  RPM  for  Unclassified 

Modelling  of  HALIFAX  Maneuvering 

Speed  (knots) 

Propeller  RPM 

5 

28.9 

10 

57.7 

15 

87.1 

20 

121.0 

25 

156.0 

30 

200.4 
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Table  9:  Rudder  Control  Properties  for  HALIFAX 


Maximum  rudder  deflection  5™dder  35  degrees 

Maximum  rudder  velocity  4™*ier  3  deg/s 

Response  natural  frequency  cog  3  rad/s 

Response  damping  Q  0.85 

Rudder-propeller  interaction  coefficient  crudder~prop  0.5 


Table  10:  Predicted  Maneuvering  Coefficients  for  HALIFAX 


Y' 

V 

-0.207 

Y' 

r 

0.062 

N' 

-0.080 

N’r 

-0.037 

y> 

Iv\v\ 

-1.006 

Y' 

Iv\r\ 

-0.140 

Y' 

■*r|r| 

0.000 

Ky 

0.000 

K\r\ 

-0.060 

NLz 

-0.200 

deflection  S^dder,  the  rudder  control  parameters  are  estimates  based  on  published 
values  for  other  ships. 

Table  10  gives  predicted  hull  maneuvering  coefficients  for  HALIFAX  determined  using 
Equations  (27)  to  (43). 

Figures  13  to  15  show  the  trajectory,  yaw  rate,  and  speed  for  HALIFAX  executing  a 
starboard  turning  circle  at  an  initial  speed  of  25  knots  and  30  degree  rudder  angle. 
Table  11  gives  associated  turning  circle  statistics. 

When  simulating  maneuvers  of  existing  ships,  hull  maneuvering  coefficients  are  often 
adjusted  to  provide  improved  agreement  between  observed  and  predicted  maneuvers. 
To  examine  the  influence  of  hull  maneuvering  coefficients  on  predicted  turning  circle 
properties,  simulations  have  been  conducted  for  which  predicted  hull  maneuvering 
coefficient  have  been  incremented  by  the  range  of  uncertainty  as  given  in  Table  1. 
Tables  12  and  13  show  the  predicted  variations  to  turning  circle  properties.  Turning 
circle  properties  are  most  responsive  to  changes  in  the  linear  yaw-yaw  maneuvering 
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Figure  13:  Turning  Circle  Trajectory  for  HALIFAX  Starboard  Turn  at  25  Knots  with 
30  Degree  Rudder 


Figure  14:  Yaw  Rate  Versus  Time  for  HALIFAX  Starboard  Turn  at  25  Knots  with 
30  Degree  Rudder 
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Figure  15:  Ship  Speed  Versus  Time  for  HALIFAX  Starboard  Turn  at  25  Knots  with 
30  Degree  Rudder 
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Table  11:  Properties  of  Predicted  HALIFAX  Starboard  Turning  Circle  at  25  Knots 
with  30  Degree  Rudder 


Transfer 

341  m 

Advance 

511  m 

Tactical  diameter 

768  m 

Final  diameter 

768  m 

Final  yaw  velocity 

1.79  m/s 

Final  speed 

23.2  knots 

coefficient  N'r.  Variation  of  the  other  yaw  force  coefficients  N'Q  and  TV',,,  also  has  a 
significant  influence  on  turning  circle  properties.  Final  speed  for  HALIFAX  varies  by 
less  than  0.5  percent  for  variation  of  each  hull  maneuvering  coefficient. 

Predicted  turning  circle  characteristics  will  also  vary  with  the  selected  value  of  the 
rudder-propeller  interaction  coefficient,  which  can  be  difficult  to  assess  accurately. 
Table  14  indicates  that  turning  circle  properties  can  vary  by  approximately  10  per¬ 
cent  when  the  rudder-propeller  interaction  coefficient  is  varied  by  0.3  relative  to  its 
reference  value  of  0.5. 
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Table  12:  Variation  of  HALIFAX  Turning  Circle  Properties  with  Increases  to  Hull 
Maneuvering  Coefficients,  25  Knots  with  30  Degree  Rudder 


Percentage  change  in 

turning  circle  property 

Maneuver 

Coef. 

Tactical 

Advance 

Transfer  Final 

Final 

Final 

coef. 

increment 

diameter 

diameter 

yaw  rate 

speed 

Y< 

V 

0.05 

0 

0 

-1 

0 

0 

0 

Y' 

-L  j. 

0.05 

1 

-1 

4 

1 

-1 

0 

N' 

1  y  V 

0.02 

7 

5 

7 

7 

-6 

0 

N’r 

0.02 

-16 

-14 

-16 

-16 

19 

0 

Y\  , 

v\v\ 

0.20 

0 

0 

-1 

0 

0 

0 

Y' 

*  v\r\ 

0.20 

0 

0 

-2 

0 

1 

0 

Y\ |  , 

r\r\ 

0.05 

0 

0 

1 

0 

0 

0 

N '  2 

vrz 

0.05 

2 

1 

2 

2 

-2 

0 

N\  , 

r\r\ 

0.02 

-6 

-4 

-6 

-6 

6 

0 

N'  2 

rvz 

0.20 

-2 

-2 

-2 

-2 

3 

0 

Table  13:  Variation  of  HALIFAX  Turning  Circle  Properties  with  Decreases  to  Hull 
Maneuvering  Coefficients,  25  Knots  with  30  Degree  Rudder 

I  Percentage  change  in  turning  circle  property 


Maneuver 

Coef. 

Tactical 

Advance 

Transfer 

Final 

Final 

Final 

coef. 

increment 

diameter 

diameter 

yaw  rate 

speed 

Y1 

V 

-0.05 

0 

0 

1 

0 

0 

0 

Y' 

-L  r 

-0.05 

-1 

1 

-3 

-1 

1 

0 

N' 

1  y  V 

-0.02 

-6 

-4 

-6 

-6 

6 

0 

N' 

r 

-0.02 

17 

14 

17 

17 

-15 

0 

Y'  , 

v\v\ 

-0.20 

0 

0 

1 

0 

0 

0 

Y' 

*  v\r\ 

-0.20 

0 

0 

1 

0 

-1 

0 

Y’ |  , 

r\r\ 

-0.05 

0 

0 

-1 

0 

0 

0 

N'  2 

vrz 

-0.05 

-2 

-1 

-2 

-2 

2 

0 

N\  , 

r\r\ 

-0.02 

5 

4 

5 

5 

-5 

0 

N ’  2 

rvz 

-0.20 

2 

2 

2 

2 

-2 

0 
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Table  14:  Variation  of  HALIFAX  Turning  Circle  Properties  with  Changes  to  Rudder- 
Propeller  Interaction  Coefficient  Crudder~prop ,  25  Knots  with  30  Degree  Rudder 

I  Percentage  change  relative  to  Crudder~prop  =  0.5 


(jr  udder— prop 

Tactical 

Advance 

Transfer 

Final 

Final 

Final 

diameter 

diameter 

yaw  rate 

speed 

0.20 

12 

9 

13 

12 

-12 

-2 

0.80 

-8 

-6 

-9 

-8 

11 

1 

30 
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13  Validation  of  Motions  in  Waves  for 
Steered  Warship  Model  on  Straight 
Course 


Previous  validation  for  ShipMo3D  [4,  24]  has  shown  very  good  agreement  with  ex¬ 
perimental  motions  for  ships  travelling  at  steady  speed  and  heading  in  waves.  The 
introduction  of  hull  maneuvering  forces  to  ShipMo3D  necessitates  new  validation  of 
predictions  in  the  frequency  and  time  domain  for  a  ship  assumed  to  have  quasi-steady 
speed  and  heading.  Furthermore,  the  new  ShipMo3D  capability  for  modelling  freely 
maneuvering  ships  requires  validation  of  predicted  motions  in  waves. 

This  section  presents  validation  of  motions  in  waves  for  the  steered  warship  model 
of  Lloyd  and  Crossland  [25],  which  was  used  for  previous  ShipMo3D  validation  in 
Reference  4.  The  validation  includes  numerical  predictions  in  the  frequency  and  time 
domains  for  a  ship  with  quasi-steady  speed  and  heading,  and  in  the  time  domain  for  a 
freely  maneuvering  ship.  All  predictions  include  the  new  ShipMo3D  hull  maneuvering 
forces.  The  predictions  for  a  freely  maneuvering  ship  also  include  forces  due  to  viscous 
resistance,  propeller  thrust,  and  rudder-propeller  interaction.  The  hull  maneuvering 
forces  are  estimated  based  on  Section  4.  For  the  predictions  assuming  quasi-steady 
speed  and  heading,  the  nonlinear  maneuvering  force  terms  are  assumed  to  be  zero. 
Table  15  gives  resistance  coefficients  for  the  steered  warship  model  estimated  using 
Schmitke  and  Murdey  [23].  The  two  propellers  are  assumed  to  have  a  diameter  of 
0.2  m  and  the  thrust  coefficient  curve  as  follows  based  on  Figure  16  of  Reference  16: 

KT(Jprop)  =  0.40  —  0.20  Jprop  —  0.15  Jprop  (81) 

Each  of  the  two  rudders  is  located  behind  a  propeller.  The  interaction  between  each 
rudder  and  its  nearby  propeller  is  modelled  using  a  propeller  interaction  coefficient 

(Jruddei  prop  yalue  Qf  Q.9. 

In  the  previous  validation  work  of  Reference  4,  extra  stiffness  and  damping  for  surge 
and  sway  were  introduced  to  ensure  stationkeeping  for  numerical  predictions  during 
low  encounter  frequencies.  In  the  present  validation,  these  extra  stiffness  and  damp¬ 
ing  terms  are  removed  to  ensure  consistency  with  experimental  conditions.  Instead, 
no  numerical  predictions  are  made  when  the  non-dimensional  encounter  frequency 
uje  \J Lj g  is  less  than  0.24. 

Figures  17  to  32  show  comparisons  between  predicted  and  experimental  motions. 
Agreement  between  predicted  and  experimental  motions  is  generally  very  good.  The 
frequency  domain  and  time  domain  predictions  assuming  quasi-steady  speed  and 
heading  give  very  similar  results.  The  predictions  for  a  freely  maneuvering  ship 
are  slightly  different  from  those  assuming  quasi-steady  speed  and  heading,  likely 
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Table  15:  Resistance  Coefficients  for  Steered  Warship  Model 


Froude  Resistance  Coefficient 

number  CR  =  FR/{l/2pU2Aw) 


<  0.150 

0.0040 

0.225 

0.0042 

0.300 

0.0057 

0.375 

0.0069 

0.450 

0.0099 

0.525 

0.0108 

because  of  the  inclusion  of  nonlinear  maneuvering  forces  and  rudder-propeller  inter¬ 
action  forces  for  the  freely  maneuvering  ship. 
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□  Model  tests 


ShipMo3D  freely  maneuvering 


—  ShipMo3D  time  domain  quasi-steady 


ShipMo3D  frequency  domain 


Figure  16:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  0  degrees,  Froude 
Number  0.28 

□  Model  tests  - ShipMo3D  time  domain  quasi-steady 

-  ShipMo3D  freely  maneuvering  .  ShipMo3D  frequency  domain 


Figure  17:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  0  degrees,  Froude 
Number  0.37 
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Yaw  RAO  1 776 1  /( ka )  Roll  RAO  |r/4|  /( ka )  Sway  RAO  1 772 1  /a 


□  Model  tests 


- ShipMo3D  time  domain  quasi-steady 

ShipMo3D  freely  maneuvering  .  ShipMo3D  frequency  domain 


Figure  18:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  30  degrees,  Froude 
Number  0.18 
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□  Model  tests 


ShipMo3D  freely  maneuvering 


—  ShipMo3D  time  domain  quasi-steady 
•  •  •  ShipMo3D  frequency  domain 


Figure  19:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  30  degrees,  Froude 
Number  0.27 
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Yaw  RAO  1 776 1  /( ka )  Roll  RAO  |r/4|  /( ka )  Sway  RAO  1 772 1  /a 


□  Model  tests 


- ShipMo3D  time  domain  quasi-steady 

ShipMo3D  freely  maneuvering  .  ShipMo3D  frequency  domain 


Wave  Frequency  loj \J L/g 


Figure  20:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  30  degrees,  Froude 
Number  0.37 
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□  Model  tests 

—  ShipMo3D  freely  maneuvering 


—  ShipMo3D  time  domain  quasi-steady 
•  •  •  ShipMo3D  frequency  domain 


Figure  21 :  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  60  degrees,  Froude 
Number  0.18 
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□  Model  tests 


- ShipMo3D  time  domain  quasi-steady 

-  ShipMo3D  freely  maneuvering  .  ShipMo3D  frequency  domain 


Wave  Frequency  ujiy/L/ g  Wave  Frequency  u>i  y/L/g 


Wave  Frequency  uii-^/L/ g  Wave  Frequency  uiy/L/g 


Figure  22:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  60  degrees,  Froude 
Number  0.27 
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Yaw  RAO  1 776 1  /( ka )  Roll  RAO  |r/4|  /( ka )  Sway  RAO  1 772 1  /a 


□  Model  tests 


- ShipMo3D  time  domain  quasi-steady 

ShipMo3D  freely  maneuvering  .  ShipMo3D  frequency  domain 


Figure  23:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  60  degrees,  Froude 
Number  0.36 
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□  Model  tests 


ShipMo3D  freely  maneuvering 


—  ShipMo3D  time  domain  quasi-steady 
•  •  •  ShipMo3D  frequency  domain 


Figure  24:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  75  degrees,  Froude 
Number  0.18 
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□  Model  tests 

—  ShipMo3D  freely  maneuvering 


—  ShipMo3D  time  domain  quasi-steady 
•  •  •  ShipMo3D  frequency  domain 


Figure  25:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  75  degrees,  Froude 
Number  0.28 
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□  Model  tests 

-  ShipMo3D  freely  maneuvering 


Wave  Frequency  uJi^/L/g 


Wave  Frequency  ujiy/L/ g 


Wave  Frequency  uii-^/L/ g 


—  ShipMo3D  time  domain  quasi-steady 
•  •  •  ShipMo3D  frequency  domain 


Wave  Frequency  uJi^L/g 


Wave  Frequency  u>i  y/L/g 


Wave  Frequency  uiy/L/g 


Figure  26:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  75  degrees,  Froude 
Number  0.36 
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□  Model  tests 


- ShipMo3D  time  domain  quasi-steady 

-  ShipMo3D  freely  maneuvering  .  ShipMo3D  frequency  domain 


Wave  Frequency  ujiy/L/ g  Wave  Frequency  uj  y/L/g 


Wave  Frequency  uii-^/L/ g  Wave  Frequency  uiy/L/g 


Figure  27:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  90  degrees,  Froude 
Number  0.18 
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□  Model  tests 


- ShipMo3D  time  domain  quasi-steady 

-  ShipMo3D  freely  maneuvering  .  ShipMo3D  frequency  domain 


Wave  Frequency  uii-^/L/ g 


Wave  Frequency  ujj\JL/ g 


Figure  28:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  90  degrees,  Froude 
Number  0.28 
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□  Model  tests 


- ShipMo3D  time  domain  quasi-steady 

ShipMo3D  freely  maneuvering  .  ShipMo3D  frequency  domain 


Wave  Frequency  ujiy/L/ g 


Figure  29:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  90  degrees,  Froude 
Number  0.36 
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□  Model  tests 


- ShipMo3D  time  domain  quasi-steady 

-  ShipMo3D  freely  maneuvering  .  ShipMo3D  frequency  domain 


Wave  Frequency  uii-^/L/ g  Wave  Frequency  uiy/L/g 


Figure  30:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  120  degrees,  Froude 
Number  0.27 
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□  Model  tests 


- ShipMo3D  time  domain  quasi-steady 

ShipMo3D  freely  maneuvering  .  ShipMo3D  frequency  domain 


S3 


o 

<1 

CD 

rT3 


Wave  Frequency  uii^jL/g 


Figure  31:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  150  degrees,  Froude 
Number  0.26 
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Yaw  RAO  1 776 1  /( ka )  Roll  RAO  |r/4|  /( ka )  Sway  RAO  1 772 1  /a 


□  Model  tests 


- ShipMo3D  time  domain  quasi-steady 

ShipMo3D  freely  maneuvering  .  ShipMo3D  frequency  domain 


Wave  Frequency  ujiy/L/ g 


Wave  Frequency  u>i  y/L/g 
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Wave  Frequency  uiy/Lj g 


Figure  32:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  180  degrees,  Froude 
Number  0.26 
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14  Recommendations 


When  computing  the  motions  of  a  freely  maneuvering  ship,  the  greatest  uncertainties 
in  forces  and  resulting  motions  are  due  to  hull  maneuvering  forces.  In  the  field  of 
maneuvering  prediction,  it  is  common  practice  to  determine  hull  maneuvering  force 
coefficients  using  model  tests.  It  is  recommended  that  future  work  examine  methods 
for  predicting  hull  maneuvering  force  coefficients.  Computational  fluid  dynamics 
(CFD)  is  likely  a  very  suitable  approach  for  estimating  hull  maneuvering  coefficients. 
Although  CFD  is  likely  too  computationally  intensive  for  direct  simulation  of  ship 
maneuvering,  prediction  of  maneuvering  coefficients  using  CFD  is  likely  feasible. 

Regression  methods  for  predicting  hull  maneuvering  coefficients,  such  as  those  pre¬ 
sented  in  Section  4,  are  much  easier  to  apply  than  CFD  predictions.  The  equations 
developed  by  Inoue  et  al.  [10]  are  based  on  data  for  commercial  ships.  Regression 
equations  based  on  data  for  naval  hull  forms  could  be  very  useful.  It  is  also  recom¬ 
mended  that  regression  approaches  be  developed  that  account  for  the  difference  in 
profiles  between  various  ship  types.  For  example,  the  shallow  draft  at  the  stern  of 
modern  naval  frigates  likely  causes  hull  maneuvering  force  coefficients  to  be  signifi¬ 
cantly  different  from  those  ships  with  constant  draft  in  the  aft  portion. 


15  Conclusions 


The  ShipMo3D  library  has  been  extended  to  model  freely  maneuvering  ships  in  calm 
water  and  in  waves.  New  force  components  have  been  added,  including  hull  maneu¬ 
vering,  hull  resistance,  propulsion,  and  rudder-propeller  interaction.  Comparisons  of 
turning  circle  predictions  with  full-scale  trials  data  for  the  tanker  Esso  Osaka  give 
encouraging  results.  Comparisons  of  predictions  with  motions  of  a  steered  warship 
model  in  waves  give  very  good  results.  Excellent  agreement  between  predictions  for 
a  freely  maneuvering  ship  and  for  a  ship  with  quasi-steady  speed  and  heading  indi¬ 
cates  that  the  extension  of  ShipMo3D  to  freely  maneuvering  ships  has  been  correctly 
implemented.  It  is  recommended  that  future  work  further  investigate  prediction  of 
hull  maneuvering  forces. 
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added  mass  matrix  for  static  foil  i 
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area  of  hull  panel  k 

added  mass  matrix  for  rudder  i 

added  mass  matrix  for  skeg  i 

hull  lateral  profile  area 

lateral  profile  area  of  section  i 

hull  wetted  surface  area 
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ship  damping  matrix 

damping  matrix  for  bilge  keel  i 

damping  matrix  for  static  foil  i 

hull  maneuvering  damping  matrix 

hull  maneuvering  damping  speed  term 
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hull  maneuvering  damping  speed  and  drift  velocity  squared  term 

hull  radiation  damping  speed  term 

hull  viscous  damping  matrix  due  to  roll 
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